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Abstract 


Effects  of  transient  rate-of-injection  (ROI)  on  jet  development  and  combustion 
were  measured  in  a  compression-ignition  direct-injection  optical  engine.  Jet  penetra¬ 
tion,  liquid  length,  jet  dispersion  angle,  liftoff  length,  and  engine  cylinder  pressure 
data  were  simultaneously  acquired  over  a  range  of  engine  operating  conditions  rel¬ 
evant  to  the  design  of  Army  engines.  Tests  were  performed  using  #2  diesel  fuel,  jet 
fuel  (JP8),  and  a  hydroprocessed  renewable  jet  fuel  (HRJ).  Ambient  thermodynamic 
conditions  and  fuel  injection  parameters  had  similar  effects  on  combustion  and  jet 
performance  in  the  optical  engine  as  in  previous  results  from  constant  volume  cham¬ 
bers.  The  higher  volatility  and  lower  density  of  the  jet  fuels  relative  to  #2  diesel  fuel 
resulted  in  shorter  liquid  lengths  for  the  jet  fuels  and  slight  differences  in  the  jet  dis¬ 
persion  and  penetration.  However,  the  lift-off  length  and  ignition  delay  for  the  JP8 
fuel  were  almost  identical  to  the  #2  diesel  fuel  with  matched  cetane  number,  indicat¬ 
ing  that  cetane  number  has  a  stronger  influence  than  changes  in  fuel  physical  prop¬ 
erties  on  the  combustion  process.  The  ramp-on  transient  from  the  opening  of  the  fuel 
injector  caused  the  scaling  of  jet  parameters  with  time  during  the  transient  ROI  to 
be  modified  from  relationships  expected  from  the  literature  for  all  of  the  fuels.  Initial 
penetration  of  the  jets  was  proportional  to  H  where  a  is  had  values  significantly  >  1 
(typically  in  the  range  of  1.5  to  2.0  dependent  upon  injection  pressure)  at  the  SOI  and 
shortly  after  in  comparison  to  the  scaling  with  time  often  assumed  in  the  literature. 
Similarities  were  seen  in  a  comparison  between  the  jet  velocity  and  the  derivative  of 
the  fuel  mass  ROI.  In  the  long-time  limit,  injections  were  seen  to  return  to  their  quasi¬ 
steady  forms  once  fuel  mass  ROI  had  become  approximately  steady,  with  the  changes 
in  jet  parameters  occurring  during  the  initial  opening  of  the  injector.  The  differences 
observed  in  jet  penetration  during  the  injector  ramp-on  transient  will  effect  the  over¬ 
all  jet  penetration  and  ignition  delay  and  liftoff  length  and  are  important  to  model, 
especially  for  injections  with  short  overall  duration  or  long  injector  transient  where 
the  transient  period  encompasses  a  significant  fraction  of  the  injection  process. 
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1  Introduction  and  Motivation 


The  compression  ignition  (Cl)  internal  combustion  engine  has  been  one  of  the  primary 
power  plants  of  the  past  100  years  due  to  its  high  efficiency  and  power-to-weight  ratio 
[1].  Most  modern  versions  use  a  fuel  delivery  system  composed  of  a  high-pressure  in¬ 
jector  mounted  to  inject  fuel  directly  into  the  combustion  chamber,  a  process  known  as 
direct  injection  (DI).  During  injection,  fuel  is  vaporized  by  and  mixed  with  the  dense  hot 
in-cylinder  gases  leading  to  autoignition  of  the  spray  This  process  is  central  to  the  perfor¬ 
mance  of  the  engine  [2],  The  start-of-combustion  is  preceded  by  a  combination  of  phys¬ 
ical  and  chemical  processes,  encompassing:  the  fluid  dynamics  of  fuel-oxidizer  mixing, 
the  thermodynamics  of  ambient  and  local  conditions  in  the  spray,  and  the  reaction  path¬ 
ways  that  ultimately  determine  the  mechanisms  of  ignition  [3].  Under  transient  injection 
conditions,  particularly  relevant  to  small  diesel  engines  of  Army  interest,  the  processes 
leading  to  ignition  and  combustion  are  directly  impacted  by  the  time-varying  rate  of  fuel 
injection.  The  current  work  has  focused  on  improving  the  understanding  of  the  transient 
period  of  the  injection  process  and  the  influence  it  has  on  the  macroscopic  behavior  of  the 
fuel  jet  up  until  ignition. 

Vaporizing  fuel  jets  have  long  been  the  subject  of  focused  study  in  the  field  of  inter¬ 
nal  combustion  engines.  The  primary  foci  of  previous  research  on  this  topic  includes:  the 
processes  governing  jet  and  spray  penetration  rates  and  distances,  flame  liftoff  lengths, 
and  ignition  location  and  delay.  The  modem  interpretation  of  the  lifted  diffusion  flame 
in  direct-injection  diesel  engines  described  by  Dec  [4]  in  1997,  and  updated  by  Pickett 
and  Siebers  [5],  involves  a  complex  coupling  of  turbulent  fluid  mechanics  with  chemical 
reactions  and  combustion,  that  remains  difficult  to  predict  through  analytical  and  compu¬ 
tational  means.  Optical  experimental  techniques  have  been  used  extensively  to  provide  a 
non-intrusive  means  for  both  gaining  a  physical  understanding  of  these  processes  and  for 
providing  quantitative  data  against  which  computational  models  can  be  vetted  or  tuned. 
A  recent  paper  by  Som  et  al.  [6]  compared  the  predictive  capabilities  of  current  state-of- 
the-art  commercial  CFD  software  packages  with  respect  to  sprays  and  combustion,  with 
the  results  of  experimental  data  taken  in  a  constant  volume  chamber  [7].  The  performance 
of  the  simulations  was  evaluated  by  their  ability  to  predict  jet  characteristics  such  as  va¬ 
por  and  liquid  penetration  of  fuel  jets,  and  combustion  parameters  including  ignition 
and  flame  stabilization.  While  the  CFD  codes  tested  were  able  to  predict  jet  and  combus¬ 
tion  behavior  during  quasi-steady  operation  relatively  well,  differences  during  transient 
periods  were  significant  [6,  7].  An  upper  limit  to  the  accuracy  of  CFD  modeling  of  fuel 
jets  is  coupled  to  the  accuracy  of  the  rate-of-injection  profile  used  for  the  fuel  boundary 
condition  [8]. 

The  bulk  of  the  research  published  to  date  has  focused  on  the  quasi-steady  period 
of  the  injection  process.  While  these  results  are  useful  for  gaining  a  fundamental  un¬ 
derstanding  of  the  combustion  and  mixing  processes  inherent  in  compression-ignition 
engines,  they  are  not  entirely  representative  of  the  actual  conditions  seen  in  real  applica¬ 
tions,  especially  in  the  case  of  the  small-bore  diesel  engine.  This  type  of  engine  often  uti¬ 
lizes  short  injection  durations  (resulting  in  ramped  rate-of-injection  profiles)  that  do  not 
correspond  to  the  conditions  under  which  many  of  the  scaling  theories  and  dependencies 
set  forth  in  the  literature  were  developed.  Additionally,  engines  operating  with  multi- 
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pie  in-cylinder  injections  per  cycle  will  similarly  have  injections  dominated  by  ramped 
rate-of-injection  profiles.  Therefore,  there  is  a  strong  need  to  gain  an  understanding  of 
transient  effects  on  jet  development  and  combustion. 

The  completed  work  addresses  this  need  through  a  fundamental  optical  study  of  tran¬ 
sient  rate-of-injection  events  in  an  optically  accessible  internal  combustion  engine.  Mea¬ 
surements  of  fuel  jet  macroscopic  behavior  were  performed  using  high-speed  imaging  of 
the  in-cylinder  fuel  injection  event  to  simultaneously  measure  liquid  penetration,  vapor¬ 
ized  fuel  penetration,  jet  spreading  angle,  and  the  liftoff  length.  Fuels  used  for  the  mea¬ 
surements  included  #2  diesel  fuel,  JP-8,  and  a  hydroprocessed  renewable  jet  fuel  (HRJ).  A 
smaller  subset  of  measurements  were  performed  with  50/50  blends  (by  volume)  of  iso- 
paraffinic  kerosene  (IPK)  and  JP-8,  and  HRJ  and  JP-8.  In-cylinder  thermodynamic  condi¬ 
tions  were  varied  at  the  start-of-injection  by  varying  the  intake  pressure  and  intake  tem¬ 
perature  to  the  engine.  Three  different  injection  pressures  were  utilized  which  provide 
three  different  rate-of-injection  ramp  profiles.  The  results  provide  measurements  over  a 
wide  range  of  conditions  relevant  to  the  design  of  Army  engines. 


2  Scientific  Background 

The  work  presented  in  this  report  focuses  on  the  development  of  transient  fuel  jets  in 
diesel  engines  and  the  coupling  between  the  transient  jet  and  the  combustion  event.  The 
following  literature  review  provides  an  overview  of  previous  work  relevant  to  the  com¬ 
pleted  project.  First,  an  overview  is  given  of  the  significant  literature  related  to  current 
understanding  of  both  quasi-steady  jets  and  recent  work  on  transient  jets,  focusing  on  liq¬ 
uid  and  vapor  penetration  behavior  and  jet  dispersion.  This  is  followed  by  a  discussion  of 
aspects  of  the  combustion  process:  jet  ignition,  flame  stabilization,  and  the  quasi-steady 
liftoff  length.  Finally,  optical  techniques  and  experimental  facilities  are  reviewed. 

The  nomenclature  used  to  describe  the  liquid  and  vapor  portions  of  the  injected  fuel 
follows  that  of  Pickett  et  al.  [5],  where  the  portion  of  the  injection  with  liquid  droplets 
present  is  referred  to  as  the  "spray,"  and  the  overall  injection  process  involving  the  spray 
as  well  as  the  region  composed  of  only  fuel  vapor  and  entrained  ambient  air  is  referred 
to  as  the  "jet." 

2.1  Quasi-Steady  and  Transient  Jet  Behavior 

Analysis  of  diesel  jets  was  already  underway  by  the  1920s  [9],  although  the  majority  of 
the  work  before  the  comprehensive  study  in  1996  by  Naber  and  Siebers  at  Sandia  Na¬ 
tional  Labs  [10]  was  performed  under  non-vaporizing  (i.e.  room  temperature)  conditions, 
focusing  on  the  development  of  experimental  correlations  for  jet  penetration  as  a  function 
of  ambient  gas  density.  Barring  a  few  notable  exceptions  [11, 12],  experimental  conditions 
for  these  analyses  were  outside  those  seen  in  internal  combustion  engines.  Also,  much 
of  the  literature  upon  which  current  understanding  of  macroscopic  jet  behavior  was  de¬ 
veloped  focuses  on  the  study  of  jets  emitted  from  injection  systems  capable  of  "top-hat" 
injection  profiles.  These  systems  display  tum-on  times  (time  for  the  injector  needle  to  fully 
open  and  reach  a  maximum  fuel  flow  rate)  on  the  order  of  50  ps  or  less.  The  goal  of  these 
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studies  in  general  was  to  study  the  quasi-steady  portion  of  the  injection  process  where  the 
rate-of-injection  is  approximately  constant,  and  the  injector  transient  is  generally  ignored. 

2.1.1  Dispersion  Angle 

The  dispersion  angle  of  a  diesel  fuel  jet  is  related  to  its  penetration  characteristics  through 
its  ambient  air  entrainment  rate.  A  larger  spray  spreading  angle  indicates  that  surround¬ 
ing  gases  are  more  quickly  entrained,  allowing  for  fasfer  vaporization  and  reduction  in 
the  total  penetration  distance.  The  dispersion  measurements  performed  by  Naber  and 
Siebers  [10]  displayed  local  spray  angles  fhaf  remained  fairly  consfanf  along  the  axis  of 
the  jet  during  the  quasi-steady  (developed)  portion  of  fhe  injection.  Based  on  simple  jet 
theory  for  a  quasi-sfeady  jef,  the  ambient  gas  entrainment  rate  is  constant  along  the  spray 
axis  and  proportional  to  three  key  variables:  injector  orifice  diamefer  (do),  fuel  injection 
velocity  {U f,  directly  proportional  to  the  square  root  of  injection/ ambient  pressure  dif¬ 
ference),  and  the  tangent  of  half  the  dispersion  angle  tan{9/2).  During  the  beginning  and 
end  of  fhe  injection  events,  an  initial  transient  period  was  noted  by  Naber  and  Siebers 
where  the  dispersion  angle  was  significantly  larger  than  the  quasi-steady  value,  an  ef¬ 
fect  the  authors  hypothesized  was  due  to  needle  movement  in  the  injector  but  was  not 
investigated  further. 

The  dispersion  angle  has  been  shown  to  be  independent  of  injecfion  pressure  in  the 
quasi-steady  regime  in  both  vaporizing  and  non-vaporizing  cases  [10].  The  spreading  an¬ 
gle  has  been  shown  by  multiple  researchers  under  non-vaporizing  cases  to  scale  with  the 
ratio  of  the  ambient  gas  density  to  the  fuel  densify  (pa/pf)  raised  fo  a  power,  defermined 
empirically  [10, 11, 12, 13, 14].  The  wide  range  of  reporfed  has  been  affribufed  fo  the  sen¬ 
sitivity  of  the  experimental  measurement  of  the  angle  to  experimental  setup  and  image 
analysis  [15,  5].  Under  vaporizing  conditions  an  overall  decrease  in  spray  dispersion  an¬ 
gles  was  observed,  although  the  effect  decreased  with  increasing  ambient  pressure.  This 
contraction  is  thought  to  be  due  to  evaporative  cooling  of  the  entrained  ambient  gases. 

Research  studying  the  effect  of  acceleration  on  unsteady  turbulent  single-phase  jets 
(notably  [16, 17])  has  shown  that  accelerating  jets  display  decreased  dispersion  angles  and 
entrainment  rates,  and  vice  versa  for  decelerating  jets.  The  change  in  entrainment  rates 
was  attributed  to  the  modified  enfrainmenf  appetife  of  the  flow,  due  fo  the  difference  in 
velocities  of  coherenf  strucfures  along  the  axial  direction  of  the  jet  [17].  This  reduction  in 
entrainment  was  modeled  by  Breidenthal  [18],  who  developed  a  self-similar  solution  for  a 
round  exponentially  accelerating  jet,  and  concluded  that  the  entrainment  rate  appetite  of 
the  large  scale  eddies  is  inhibited  if  they  are  followed  by  upstream  neighboring  structures 
of  greafer  momenfum.  Other  CFD  studies  were  carried  out  corroborating  these  findings, 
including  the  unsteady  Reynolds  averaged  Navier-Stokes  simulations  of  Abani  and  Reitz 
[19]  and  direct  numerical  simulations  (DNS)  by  Silva  et  al.  [20].  The  dispersion  angle 
was  measured  during  the  accelerating  portion  of  a  non-vaporizing  jef  by  Payri  et.  al.  [21], 
where  the  angle  decreased  from  ifs  original  value  until  the  jet  reached  a  quasi-steady  state, 
similar  to  the  results  seen  from  Naber  and  Siebers  [10].  Little  to  no  experimental  data  on 
this  parameter  currently  exists  for  fransienf  jefs  under  diesel  engine-like  conditions. 
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2.1.2  Jet  Penetration 


The  propagation  of  the  head  of  the  fuel  jet  during  the  injection  process  is  defined  as  the 
jet  penetration.  Naber  and  Siebers  [?]  were  able  to  collapse  all  of  their  spray  penetration 
data  for  all  conditions  using  non-dimensional  time  and  distance  parameters  derived  from 
a  simple  jet  model  based  on  the  earlier  work  of  Wakuri  et  al.  [22],  The  jet  model  assumes  a 
uniform  steady  velocity  profile  at  any  downstream  position  x,  and  a  constant  dispersion 
angle  a.  A  no-slip  condition  was  assumed  between  the  fuel  in  the  spray  and  the  entrained 
air,  and  it  was  also  assumed  that  the  injection  velocity  was  constant  and  neglects  the  initial 
tum-on  transient.  The  analysis  is  technically  only  valid  for  incompressible  and  isothermal 
jets,  although  the  authors  successfully  applied  it  to  two-phase  jets  under  vaporizing  con¬ 
ditions  as  well. 

Naber  and  Siebers  ignored  the  initial  rate-of-injection  (ROI)  ramp-up  portion  by  ex¬ 
trapolating  the  "steady"  ROI  portion  of  the  curve  to  zero  penetration  and  defining  that 
point  as  time  zero.  Defining  t=0  in  this  way  gave  better  conformation  of  the  experimental 
data  to  the  model,  although  it  removes  the  ability  of  the  model  to  predict  the  effect  of 
initial  transients.  The  final  penetration  correlation  is  provided  in  Equations  1  and  2  in  the 
short  and  long  time  limits,  respectively: 


where  S  is  the  penetration  distance  of  the  jet,  dependent  on  the  nozzle  velocity  and  area 
coefficients  Cy  and  Ca,  pressure  drop  across  the  orifice  AP,  ambient  and  fuel  densities  pa 
and  p f,  constant  a,  orifice  diameter  do,  and  time  t.  These  correlations  show  the  same  de¬ 
pendencies  on  pressure  drop  and  time  as  previous  work  by  Hiroyasu  and  Arai  [11],  Dent 
[23],  and  others.  Also,  the  long-time  scaling  dependencies  are  the  same  as  those  seen  for 
steady-state  turbulent  gas  jets  [24].  Naber  and  Siebers  also  noted  that  when  substituting 
typical  orifice  parameters  into  their  correlation,  the  initial  scaling  factor  is  also  similar  to 
previous  work  (2.95  for  Hiroyasu  and  Arai,  3.07  for  Dent,  and  2.9  for  Naber  and  Siebers) 
[10].  The  transition  period  between  these  short  and  long-time  behaviors  is  referred  to  in 
Hiroyasu  and  Arai  as  the  spray  breakup  time,  given  in  Equation  3, 

break  =  28.65  (3) 

V  PaAp 

and  occurs  as  the  jet  shifts  from  being  dominated  by  the  liquid  to  the  entrained  gas. 

These  correlations  identify  a  linear  dependence  of  penetration  on  time  for  short  times 
and  a  square  root  dependence  for  long  times.  Recent  work  [8]  suggests,  however,  that 
the  simple  assumption  of  linear  penetration  at  injector  turn  on  is  not  necessarily  a  good 
one  and  that  the  behavior  may  be  more  complex  even  for  injectors  with  near  top-hat 
profiles.  Top-hat  injection  profiles  assume  that  the  liquid  spray  velocity  at  the  orifice  is 
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immediately  at  its  maximum  value  as  soon  as  the  injection  event  occurs,  negating  the 
initial  ramp-up  period. 

Relatively  little  data  have  been  gathered  on  vaporizing  transient  jets  under  diesel  en¬ 
gine  conditions.  The  scaling  of  spray  penetration  for  non-vaporizing  jefs  during  fhe  open¬ 
ing  of  fhe  injecfor  was  sfudied  by  Kosfas  ef  al.  [25]  for  ambienf  pressures  up  fo  5  MPa 
using  a  single-hole  common-rail  diesel  injecfor  wifh  high-speed  imaging.  The  aufhors 
found  fhaf  for  the  initial  period  during  the  opening  of  the  injector,  the  injection  velocity 
appeared  to  scale  with  time  as  U  (x  y/i,  while  the  penetration  distance  scaled  with  time 
as  S'  oc  Recent  work  by  Pickett  et  al.  [8]  has  provided  initial  insight  into  the  effect  of 
transienf  rafe-of-injecfion  on  overall  jef  developmenf.  Analysis  of  mass-  and  momenfum- 
based  ROI  dafa  revealed  fhaf  fhe  initial  fransienf  opening  time  depends  on  fhe  fofal  flow 
area  in  fhe  injecfor  and  fhaf  significanf  flucfuafions  in  the  ROI  occur  even  after  the  injector 
appears  to  have  fully  opened. 

Efforfs  have  been  made  fo  model  fransienf  jef  behavior  with  a  limited  computational 
cost  using  one-dimensional  models  of  non-reacfing,  non-vaporizing  jefs  [26,  27].  This 
model  was  recenfly  updafed  by  Pickeff  ef.  al.  [5]  fo  allow  for  variable  dispersion  angles. 
Employing  fhe  experimenfal  ROI  for  fhe  specific  injecfor  used  in  fhe  measuremenfs  and 
allowing  for  a  variable  jef  spreading  angle  greafly  enables  close  mafching  of  experimenfal 
dafa  [5]. 

2.1.3  Liquid  Length 

The  liquid  length  of  a  jef  refers  fo  fhe  disfance  from  fhe  injecfor  orifice  fo  fhe  poinf  where 
liquid  droplefs  are  no  longer  defecfable.  Under  high-femperafure  vaporizing  conditions, 
such  as  those  found  in  infernal  combustion  engines,  the  jet  begins  as  a  spray  where  liquid 
penetration  is  equal  to  the  overall  jet  penetration.  As  hot  ambient  gases  are  mixed  into  the 
spray,  a  point  in  the  jet  is  reached  where  sufficient  energy  for  all  droplefs  fo  evaporafe  is 
entrained  and  the  liquid  phase  completely  disappears.  Eurther  penetration  at  this  point 
is  entirely  vapor,  and  the  liquid  length  fluctuates  due  to  turbulence  about  a  quasi-steady 
position,  if  ambienf  condifions  remain  consfanf  [28]. 

The  effecfs  of  ambienf  femperafure  and  densify  on  fhe  liquid  length  of  vaporizing  jefs 
were  explored  by  Espey  and  Dec  [29].  Higher  femperafures  provide  more  fhermal  energy 
for  vaporization,  decreasing  liquid  lengfh  and  reducing  chemical  timescales  for  ignition. 
Higher  in-cylinder  densifies  increase  ambienf  air  enfrainment  info  fhe  jef.  This  leads  fo  an 
enfrainmenf -induced  increase  in  fhe  vaporizafion  rafe. 

Liquid  lengfhs  decrease  linearly  with  injector  orifice  diamefer,  approaching  zero  as 
orifice  diamefer  is  decreased  towards  zero  [30,  28].  A  fuel  jef 's  liquid  lengfh  is  essentially 
independenf  of  injection  pressure[10,  31]  since  fhe  increased  rafe  of  injected  fuel  mass  is 
balanced  by  a  corresponding  increase  ambient  gase  entrainment  rate.  The  vaporization  of 
fuel  sprays  under  diesel  engine  conditions  are  dominated  by  entrainment  through  turbu¬ 
lent  mixing,  rather  than  "local  interphase  transport  mechanics"  [28].  This  does  not  mean 
that  these  other  processes  are  not  important,  only  that  they  are  not  rate-limiting  under 
typical  diesel  engine  injection  conditions. 

Liquid  length  is  non-linearly  dependent  on  ambient  density  and  temperature.  It  also 
is  affected  by  fuel  properties,  namely  fuel  volafilify.  The  effecfs  of  fuel  fype  on  liquid 
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length  have  been  documented  in  multiple  works,  including  [28,  32,  33].  Higgins  [33]  also 
developed  an  engineering  correlation  from  empirical  data  to  predict  liquid  length  for  a 
wide  range  of  fuels,  given  a  set  of  fuel  properties.  These  correlations  are  based  on  results 
from  quasi-steady  injection  measurements. 

Fisher  and  Mueller  [34,  35,  36]  recently  examined  effects  of  transient  conditions  on 
fuel  sprays  in  an  optical  engine.  The  experimental  results  (with  the  start  of  injection  taken 
over  a  range  of  engine  timings)  matched  those  from  Sieber's  mixing  limited  vaporiza¬ 
tion  model  [30],  as  long  as  the  instantaneous  thermodynamic  condition  were  used  in  the 
model.  Increased  temperature  due  to  combustion  heat  release  reduced  the  mean  liquid 
length. 

2.2  Liftoff  Length  and  Ignition 

The  mixing  and  vaporization  processes  that  dominate  the  formation  of  jets  as  described  in 
the  previous  section  also  serve  as  determining  factors  for  ignition  and  flame  stabilization 
in  the  engine.  The  conceptual  diesel  jet  combustion  model,  as  developed  by  Dec  et  al. 
[4,  37]  and  more  recently  updated  by  Pickett  and  Siebers  [38],  describes  the  processes 
for  mixing-controlled  diesel  combustion  where  a  rich  partially-premixed  air-fuel  mixture 
begins  to  combust  at  a  downstream  distance  referred  to  as  the  liftoff  length,  creating  a  rich 
high-temperature  region  (suitable  for  soot  formation)  surrounded  by  a  turbulent  non- 
premixed  flame.  After  the  initial  autoignition  of  the  fuel  jet,  a  quasi-steady  liftoff  length  is 
established  that  is  influenced  by  changes  in  the  in-cylinder  conditions.  The  exact  physics 
behind  the  establishment  of  this  location  are  still  not  fully  known,  and  the  focus  of  much 
of  the  work  performed  on  combusting  jets  has  been  devoted  to  discerning  whether  this  is 
an  ignition-based  or  flame-propagation-driven  process. 

In  hydrocarbon-oxygen  chemical  reactions,  significant  OH  chemiluminescence  occurs 
during  high-temperature,  near-stoichiometric  combustion  [39],  and  the  beginning  of  this 
region  marks  the  liftoff  length  of  the  flame  [40].  This  emission  comes  from  the  excited 
hydroxyl  radical,  OH*,  formed  during  the  exothermic  reaction  [39]: 

CH  +  02^00  +  OH*.  (4) 

The  radical  quickly  returns  to  its  ground  state  through  a  combination  of  collisional  quench¬ 
ing  and  emission  of  chemiluminescence,  the  strongest  band  of  which  lies  in  a  narrow 
region  around  310  nm  [39]. 

Siebers  and  Higgins  [41]  used  emision  from  OH*  to  measure  quasi-steady  liftoff  lengths 
in  a  CVC  over  a  large  range  of  ambient  thermodynamic  and  injection  conditions  inclu¬ 
sive  of  those  seen  in  modern  diesel  engines  and  determined  power-law  dependencies  for 
liftoff  on  each  parameter.  The  scaling  relations  were  used  to  develop  a  power  law  corre¬ 
lation  [42]: 

H  =  (5) 

where  H  [mm]  is  the  liftoff  length,  O  is  a  proportionality  constant,  [K]  is  the  ambient 
gas  temperature,  pa  [kg/m^]  is  the  ambient  gas  density,  do  [mm]  is  the  injector  orifice  di¬ 
ameter,  U  [m/ s]  is  the  fuel  injection  velocity  at  the  orifice,  and  Zst  is  the  stoichiometric 
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mixture  fraction.  Liftoff  length  exhibits  a  strong  non-linear  dependence  on  thermody¬ 
namic  conditions,  an  effect  that  the  authors  postulated  is  due  in  part  to  their  influence 
on  the  laminar  flame  speed  and  thermal  diffusivity.  Also,  unlike  liquid  length,  the  liftoff 
length  is  linearly  dependent  upon  injection  velocity,  a  behavior  seen  in  combusting  gas 
jets  [40]. 

Fuel  autoignition  chemistry  also  directly  effects  the  liftoff  length.  Fuels  with  the  longer 
ignition  delays  generally  show  longer  liftoff  lengths  [42].  This  potentially  indicates  that 
ignition  is  indeed  a  factor  in  the  stabilization  of  the  lifted  flame.  However,  there  are  cases 
where  the  liftoff  length  ordering  does  not  follow  ordering  based  on  ignition  delay,  indicat¬ 
ing  that  the  correlation  between  ignition  delay  and  liftoff  length  is  not  always  one-to-one 
[42].  Additional  liftoff  fuels  studies  include  [43, 44, 45]. 

It  is  specifically  mentioned  by  Pickett  et  al.  [46]  that  the  correlation  given  in  Equation 
5  likely  would  not  hold  for  transient  jets,  due  to  the  coupling  of  the  spatial  location  of  the 
initial  flame  kernels  with  the  established  quasi-steady  location.  CFD  simulations  were 
performed  by  Juneja  et  al.  [47]  exploring  the  effect  of  the  injection  rate-shape  on  diesel 
combustion  where  jets  with  top-hat,  accelerating,  and  decelerating  ROIs  were  tested.  The 
modeling  was  performed  using  detailed  chemistry  in  a  simulated  CVC  that  was  able  to 
produce  lifted-flame  behavior.  The  top-hat  ROI  profile  simulations  agreed  with  the  results 
from  Pickett  et  al.  [42],  but  significant  differences  were  seen  in  the  cases  with  accelerat¬ 
ing  and  decelerating  ROI  profiles.  For  the  decelerating  rate  shape,  the  liftoff  length  was 
shown  to  decrease  over  time,  an  expected  result  as  the  injection  velocity  was  decreasing. 
However,  for  the  accelerating  case,  the  initial  liftoff  length  was  short  due  to  the  initial  low 
injection  velocity  and  was  maintained  at  this  short  length  over  the  entire  injection  (a  4  ms 
event).  It  was  postulated  by  the  authors  that  this  was  due  to  hot  combustion  gases  being 
built  up  around  the  initial  ignition  zone  that  then  continued  to  convect  into  the  mixing 
jet,  such  that  even  as  the  injection  velocity  increased  by  a  factor  of  two,  the  liftoff  length 
did  not  change  appreciably.  Regardless  of  the  reasons,  it  was  evident  that  transients  in  the 
ROI  altered  the  liftoff  length  from  its  normal  quasi-steady  value. 

2.3  Optics  and  Image  Analysis 

The  following  section  provides  a  brief  review  of  optical  techniques  relevant  to  sprays  and 
combusting  jets,  focused  predominantly  on  techniques  that  provide  path  integrated  mea¬ 
surements  that  can  be  performed  at  very  high  repetition  rates  (limited  mainly  by  camera 
technology)  allowing  for  the  temporal  dynamics  of  diesel  jets  to  be  resolved.  These  tech¬ 
niques  include  schlieren  and  shadowgraphy  for  vapor  measurements,  Mie  scattering  for 
liquid  penetration  measurements,  and  chemiluminescence  for  liftoff  length  measurement. 

Schlieren  and  shadowgraph  techniques  rely  upon  the  change  in  refractive  index  asso¬ 
ciated  with  density  gradients  in  a  gases  to  bend  light  resulting  in  light  and  dark  regions  in 
an  image  [48].  A  linear  relationship  exists  between  the  density  of  a  gas  and  its  refractive 
index,  given  by: 

n  —  I  =  kp  (6) 

where  n  is  the  refractive  index  of  the  gas,  k  is  the  Gladstone-Dale  constant,  and  p  is  the 
density  of  the  gas. 
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In  shadowgraphy,  light  from  a  source  is  refracted  as  it  moves  through  a  gas  with  a 
density  gradient.  Without  the  "shadowgraphy  object"  S  in  the  way  of  the  light  source,  the 
screen  would  be  uniformly  lit,  but  the  gradients  in  the  gas  cause  shadows  to  be  cast  in 
the  far  field,  produced  as  the  direction  of  the  rays  are  modified  by  the  test  region  [48].  In 
a  simple  schlieren  setup,  a  collimated  light  source  is  passed  through  a  test  region  and  the 
light  is  refocused  through  a  knife  edge  before  being  imaged  onto  a  surface.  The  knife  edge 
blocks  light  rays  that  were  refracted  by  the  density  gradients,  causing  the  re-expanded 
image  to  have  dark  regions  where  density  gradients  exist  in  the  gas  [48]. 

Chemiluminescence  and  incandescence  that  occur  as  a  byproduct  of  combustion  with¬ 
out  external  laser  excitation  are  also  useful  sources  for  optical  imaging.  Chemilumines¬ 
cence  is  radiation  emitted  from  an  excited  species  formed  in  a  chemical  reaction  [49], 
while  incandescence  is  thermal  radiation  from  heated  particles.  Chemiluminescence  for 
a  specific  transition  of  a  specific  species  occurs  in  a  specified  wavelength  band,  and  can 
therefore  be  targeted  and  imaged  using  a  bandpass  filter  on  the  camera  system.  Incandes¬ 
cence  on  the  other  hand  is  broadband,  exhibiting  a  Planck  distribution  as  modified  by  the 
spectral  emissivity  of  the  emitting  substance. 


3  Experimental  Method 

The  primary  experimental  apparatus  used  for  the  measurements  is  composed  of  a  cus¬ 
tomized  single-cylinder  compression-ignition  optical  engine  with  engine  speed  controlled 
by  a  DC  dynamometer.  The  engine  is  paired  with  a  high-speed  data  acquisition  system 
for  engine  cylinder-pressure  measurements,  a  high-pressure  common-rail  fuel  injection 
system,  and  an  optical  system  designed  to  resolve  the  jet  and  combustion  characteristics 
during  and  after  the  ROI  ramp.  The  following  is  a  description  of  the  systems  used  exper¬ 
iments,  as  well  as,  a  discussion  of  the  post-processing  algorithm  used  in  analyzing  the 
optical  data. 

3.1  Optical  Engine 

The  engine  block  for  the  optical  engine  is  a  single  cylinder  engine  block  with  balance 
shafts.  A  custom  optically  accessible  cylinder  head  and  piston  were  used  for  the  current 
work.  Schematic  cross  sections  of  the  engine  components  used  in  the  experiments  are 
shown  in  Figure  1. 

General  mechanical  specifications  of  the  engine  are  summarized  in  Table  1.  The  cylin¬ 
der  head  is  an  oil-less  design  that  provides  direct  access  to  valvetrain  components,  has 
centrally  located  fuel  injector  access,  and  features  overhead  optical  access  that  covers  ap¬ 
proximately  25%  of  the  combustion  chamber  and  extends  to  the  wall  to  allow  for  jet-wall 
boundary  analysis.  For  the  current  work,  a  stainless  steel  mirror  optically  polished  to 
<20  A  was  installed  in  head  window  location  to  allow  for  shadowgraphy  to  be  simul¬ 
taneously  performed  with  Mie  scattering  measurements  (process  detailed  below).  The 
intake  valve  is  shrouded,  similar  to  the  design  in  [50],  which  allows  for  control  over  the 
in-cylinder  swirl  ratio.  The  intake  system  for  the  engine  provides  control  of  in-cylinder 
density  through  control  of  the  mass  flow  rate  with  a  system  of  choked  flow  orifices  of 
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Figure  1:  Schematic  cross-section  of  the  optical  engine  used  in  work-to-date. 


Table  1:  Engine  Specifications 


Parameter 

Value 

Engine  type 

Single-cylinder  optical 

Bore 

82  [mm] 

Stroke 

76.2  [mm] 

Piston  optical  access  diameter 

53.3  [mm] 

Geometric  compression  ratio 

14.03:1 

Connecting  rod  length 

144.8  [mm] 

Swirl  ratio  range 

0-5.5 

Number  of  exhaust  valves/intake  valves 

1/1 

Cycle 

4-stroke 

EVO/IVO/EVC/IVC  (degrees  after  compression  TDC) 

164/346/374/556  [CAD] 

various  sizes.  Cylinder  pressure  is  measured  using  a  piezoelectric  transducer  (Kistler, 
605BB60)  sampled  at  0.25  CAD  resolution. 

This  engine  utilizes  the  Bow  ditch  extended  piston  design  [51],  which  allows  for  optical 
access  through  a  window  in  the  piston.  The  piston  crown  holds  a  29.0  mm  thick  UV 
grade  fused  silica  window  fhaf  provides  a  fop-dead  cenfer  (TDC)  field-of-view  of  53.3 
mm.  The  pisfon  bowl  is  a  righf  cylinder,  whose  size  was  designed  to  allow  for  maximum 
length  of  jet  penetration,  while  maintaining  a  sufficient  compression  ratio  for  compression 
ignition.  Due  to  the  proximity  of  the  jet  to  the  optical  window  surface,  the  engine  had  to 
be  disassembled  and  cleaned  after  every  three  injections  due  to  soot  buildup  on  the  piston 
window. 

3.2  Fuel  Injection  System 

The  fuel  injection  sysfem  is  composed  of  a  second  generation  Bosch  common-rail  setup 
capable  of  sfable  injection  pressures  from  25-160  MPa.  A  Bosch  CRIN  2  injeefor  outfitted 
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with  a  custom,  6-hole  injector  tip  from  Stanadyne  Inc.  was  used.  With  traditional  common 
rail  injectors,  total  flow  area  is  one  of  the  two  methods  of  controlling  the  transient  period 
of  the  injection,  the  other  being  the  injection  pressure.  Increasing  the  total  nozzle  flow  area 
or  decreasing  the  injection  pressure  increases  the  ramp-up  period.  This  type  of  injector 
operates  by  using  the  fuel  injection  pressure  to  overcome  the  force  of  a  spring  holding 
the  tip  needle  closed.  Higher  fuel  pressures  open  this  spring  faster,  reducing  transient 
period.  To  maximize  the  temporal  resolution  of  the  data  over  the  transient  period,  a  six- 
hole  injector  was  used  (as  compared  to  similar  tips  that  were  available  with  only  two  or 
four  tips),  and  injection  pressure  was  a  swept  variable  in  the  experiments  to  vary  transient 
time. 

The  jet  exits  the  injector  at  an  angle  of  13°  out  of  the  imaging  plane,  giving  a  maximum 
jet  penetration  distance  along  the  spray  axis  of  29.4  mm  before  impinging  upon  the  piston 
bowl.  The  hydro-ground  orifices  have  diameters  of  110  /im  with  k-factor  =  1.5,  represen¬ 
tative  of  injectors  used  in  small  compression-ignition  engines.  Injections  were  controlled 
electronically,  with  injection  timing  determined  from  the  optical  encoder  coupled  to  the 
engine  crankshaft.  The  flow  properties  of  the  injector  and  the  various  tips  have  been  char¬ 
acterized  using  a  Bosch  rate-of-injection  meter  [52]  outfitted  with  a  piezoelectric  pressure 
transducer  for  rate-shape  measurements  [53].  The  results  comparing  ROI  of  the  maximum 
and  minimum  injection  pressures  during  the  initial  transient  portion  of  the  rate  shape  are 
shown  in  Figure  2. 


Figure  2:  Rate-of-injection  during  the  transient  turn-on  period  for  maximum  and  mini¬ 
mum  injection  pressures. 

The  dependence  of  the  rate  of  increase  of  injection  velocity  on  the  injection  pressure  is 
visible  in  these  plots.  The  range  of  injecfion  pressures  examined  allows  for  durafions  of 
the  initial  injection  transient  of  up  to  500  /is. 
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Table  2:  Engine  Specifications 


Property 

Units 

#2  diesel 

JP-8  POSF  6169 

HRJ 

H/ C  ratio  (mole  based) 

1.77 

1.94 

2.15 

Net  Heat  of  Combustion 

MJ/kg 

42.1 

43.4 

44.1 

Cetane  Number 

46.1 

44.8 

57.4 

Density  (15°C) 

g/ cm^ 

0.865 

0.798 

0.739 

Viscosity  (40°C) 

cSt 

2.71 

1.30 

0.75 

3.3  Fuels  and  Engine  Operating  Conditions 

Three  primary  fuels  were  used  for  the  reported  work.  The  first  was  an  ultra-low  sulfur 
#2  diesel  fuel,  chosen  due  to  its  prevalence  in  the  literature  and  its  ubiquity  in  modem 
compression  ignition  engines.  A  JP-8  fuel  and  a  hydro-processed  renewable  jet  fuel  were 
chosen  as  the  remaining  test  fuels  due  to  their  relevance  for  Army  applications  and  to 
enable  the  influence  of  fuel  properties  to  be  investigated.  The  jet  fuels  have  lower  viscosi¬ 
ties  and  densities  than  the  diesel,  as  well  as  higher  volatilities,  as  seen  in  Figure  3.  The 
physical  properties  of  the  fuels  are  provided  in  Table  2. 


Figure  3:  Distillation  curves  for  #2  diesel  fuel,  JP-8  POSF  6169  and  FIRJ. 

Each  fuel  was  measured  at  the  operating  conditions  given  in  Table  3.  Measurements 
were  performed  at  two  in-cylinder  densities  and  three  in-cylinder  temperatures.  Three 
different  injection  pressures  were  tested  at  each  density/ temperature  condition.  This  sweep 
of  injection  pressure  provided  a  range  of  transient  periods  during  the  opening  of  the  in¬ 
jector  nozzle.  Tower  injection  pressures  result  in  longer  ramp-up  times,  while  increasing 
injection  pressure  results  in  ROI  profiles  approaching  the  square  wave  profile  used  in 
much  of  the  literature.  The  in-cylinder  swirl  ratio  was  set  to  zero  for  all  cases. 
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Table  3:  Experimental  Test  Program 


Property 

Units 

Value 

In-Cylinder  Density  @  TDC 

kg/m^ 

15,23 

Fuel  Injection  Pressure 

MPa 

30,  52.5, 75 

Injection  Duration 

US 

1500 

SOI  Timing 

[CAD] 

-5 

Engine  Speed 

RPM 

1200 

Intake  Air  Temperature 

°C 

115, 125, 135 

3.4  Optics  and  Imaging 

The  high-speed  optical  imaging  system  used  for  the  work  to  date  is  capable  of  imaging 
the  entire  injection  and  combustion  event  while  maintaining  sufficiently  high  temporal 
resolution  to  resolve  much  of  the  initial  fuel  injection  transient.  The  system  was  designed 
to  provide  simultaneous  imaging  of  the  vapor  phase  of  the  jet,  the  liquid  phase,  and  the 
region  of  high-temperature  combustion  in  the  downstream  portion  of  the  jet  in  an  internal 
combustion  engine.  In  contrast,  the  majority  of  the  relevant  literature  either  images  in  a 
constant  volume  chamber  or  captures  these  processes  separately,  not  simultaneously.  The 
system  layout  is  shown  in  two  views  in  Figure  4. 
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Figure  4:  High-speed  imaging  system  used  for  optical  measurements  of  in-cylinder  pro¬ 
cesses.  Top  view  on  the  right,  side  view  on  the  left. 


Direction  of  light  into  and  out  of  the  engine  is  controlled  by  a  pair  of  beamsplitters. 
A  longpass  dichroic  beamsplitter  (Semrock,  FF347-Di01-50.4x71.2)  reflects  the  ultravio¬ 
let  hydroxyl  radical  (OH)  chemiluminescence,  an  indicator  of  high  temperature  combus¬ 
tion  and  therefore  the  liftoff  length,  and  transmits  light  in  the  visible  spectrum  including 
light  from  the  LED  used  for  shadowgraphy  and  Mie  scattering  illumination.  A  second, 
broadband  beamsplitter  (Edmund  Optics,  NT48-900)  reflects  30%  of  the  incident  light 
and  transmits  70%  as  seen  on  the  left  in  Figure  4.  The  combustion  chamber  is  illuminated 
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by  a  high  power  forced  air-cooled  LED  (with  maximum  light  output  of  7.2  W)  emitting 
light  at  a  center  wavelength  at  462  nm  (Luminus,  CBT-120)  that  provides  the  lighting  for 
jet  visualization.  A  set  of  collimating  optics  are  used  to  introduce  the  light  emitted  from 
the  LED  into  the  engine.  The  LED  light  that  is  immediately  lost  as  it  is  reflected  by  the 
broadband  visible  beamsplitter  is  caught  by  an  absorbing  surface  to  prevent  uimecessary 
reflections.  Although  the  LED  must  be  run  at  a  higher  power  than  would  typically  be 
required  due  to  this  initial  70%  loss,  it  is  not  a  limiting  factor  due  to  the  high  power  ca¬ 
pabilities  of  the  LED  system.  The  LED  itself  is  a  l"xl"  square  circuit  board  with  a  small 
cooling  system  and  remotely  located  power  supply.  A  non-coherent  light  source  such  as 
an  LED  does  not  exhibit  speckle  [54]  seen  in  more  traditional  laser-based  setups. 

The  camera  (Vision  Research,  Phantom  v311)  used  for  jet  visualization  is  operated  at 
a  frame  rate  of  120.1  kHz  with  an  exposure  duration  of  2.0  ps.  The  camera  is  outfitted 
with  a  E/1.4  85  mm  focal  length  Nikkor  lens.  This  frame  rate  provides  for  the  capture 
of  60  frames  within  the  initial  500  ps  transient,  allowing  for  resolution  of  the  jet  propa¬ 
gation  and  development  for  analysis  while  still  maintaining  a  spatial  resolution  of  256 
X  80  square  pixels,  at  a  pixel  side  size  of  125  pm  in  the  field-of-view.  The  low  exposure 
duration  reduces  motion  blur  by  limiting  movement  of  the  jet  during  a  single  frame,  ef¬ 
fectively  freezing  the  flow  motion  over  most  of  the  field-of-view.  A  60  nm  bandpass  filter 
centered  at  460  nm  is  used  to  help  reject  much  of  the  broadband  soot  incandescence  and 
combustion  luminosity  while  transmitting  the  LED  light  scattered  back  from  the  engine. 

A  second  high-speed  camera  (Vision  Research,  Phantom  v7.1)  outfitted  with  a  lens- 
coupled  high-speed  intensifier  (Night-Vision  Systems,  S20  photocathode  with  P46  phos¬ 
phor)  is  used  tor  OH  chemiluminescence  imaging  and  liftoff  length  measurements.  OH 
chemiluminescence  centered  in  the  UV  at  310  nm  is  imaged  onto  the  intensifier  through  a 
combination  of  a  280  nm  longpass  filter  (Schott  Glass,  WG280)  and  a  80  nm  PWHM  band¬ 
pass  filter  centered  at  300  nm  (Semrock,  PPOl-300/80-25)  fitted  to  a  UV  P/4.5  105  mm 
focal  length  Nikkor  lens.  The  camera  is  operated  at  a  frame  rate  of  40  kHz  and  an  expo¬ 
sure  duration  of  20  ps.  The  longer  exposure  duration  and  slower  frequency  were  required 
to  obtain  enough  OH  chemiluminescence  signal  for  images  of  sufficient  quality.  Also,  it  is 
less  necessary  for  extremely  high  (>100  kHz)  frame  rates  in  the  region  where  this  camera 
sees  signals  due  to  the  downstream  location  of  the  lifted  flame  and  concomitant  lower 
local  velocities.  This  method  of  imaging  the  location  of  high-temperature  combustion  via 
the  chemiluminescence  of  the  excited  OH  radical  is  well  documented  for  use  in  optical 
engines  [44, 55].  The  vapor  and  liquid  techniques  used  here  are  less  conventional,  and  are 
explained  briefly. 

3.4.1  Vapor  Imaging:  Shadowgraphy 

Imaging  transparent  gases  such  as  fuel  vapors  requires  special  techniques  such  as  Schlieren 
or  shadowgraphy.  As  described  in  the  scientific  background,  the  index  of  refraction  of  a 
gas  is  a  function  of  its  chemical  makeup,  density,  and  temperature.  Rays  of  light  passing 
through  regions  that  have  gradients  of  these  properties  are  steered,  and  can  be  focused 
into  an  image  that  has  light  and  dark  regions  corresponding  to  the  location  of  those  gra¬ 
dients.  The  polished  metal  insert  in  the  cylinder  head  allows  for  double-pass  shadowgra¬ 
phy  to  be  used  to  detect  vapor  penetration.  This  surface  promotes  specular  reflection  of 
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the  incoming  LED  light.  As  the  incident  and  reflected  light  propagates  through  the  fuel 
jet,  gradients  in  density  refract  the  light  and  change  the  amount  of  light  returning  to  the 
camera,  providing  a  clear  indicator  for  the  jet  vapor  penetration.  A  diagram  illustrating 
the  process  via  ray  tracing  is  shown  in  Figure  5.  This  is  a  graphic  representing  the  actual 
ray- tracing  setup  built  in  a  CAD  program  that  determines  the  quantitative  effects  of  the 
presence  of  density  and  temperature  gradients  due  to  the  refraction  of  the  incoming  LED 
light  from  the  fuel  vapor. 
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Figure  5:  A  ray-tracing  diagram  of  the  paths  followed  by  the  light  imaged  for  shadowg- 
raphy  in  vapor  measurements. 

The  collimated  LED  light  is  introduced  at  angle  of  incidence  of  (1-5°)  relative  to  the 
the  head  mirror,  such  that  when  no  density  gradients  exist  in-cylinder,  the  reflected  light 
returns  on  a  path  that  does  not  reach  the  camera  sensor,  and  the  image  of  the  mirror 
appears  dark.  When  the  fuel  jet  passes  in  front  of  the  mirror,  providing  gradients  in  den¬ 
sity  and  temperature,  the  light  rays  are  refracted  to  an  angle  closer  to  perpendicular  with 
the  mirror  surface,  and  return  along  a  path  to  the  camera  sensor.  The  term  "double-pass 
shadow graphy"  comes  from  the  fact  that  the  beam  is  initially  steered  as  it  passes  through 
the  jet  as  it  initially  comes  into  the  combustion  chamber,  and  again  after  reflecting  off  of 
the  mirrored  surface.  The  error  in  position  introduced  from  the  short  distance  (<  5  mm) 
traveled  between  the  jet  and  the  mirror  surface  is  negligible. 

A  similar  setup  for  imaging  liquid  and  vapor  penetration  was  used  in  [56],  although  in 
that  work  the  hquid  and  vapor  diagnostics  had  to  be  imaged  separately  and  chemilumi¬ 
nescence  was  not  captured.  Measurements  involving  fuel  vapor  usually  involve  Schlieren 
techniques  which  require  direct  optical  access  on  a  linear  path  into  and  then  out  of  the 
apparatus,  access  which  is  possible  with  a  CVC  [5]  but  typically  impossible  in  an  engine 
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setting. 


3.4.2  Liquid  Imaging:  Extinction 

The  liquid  portion  of  the  spray  is  composed  of  a  dense  region  of  very  small  droplefs  with 
size  on  the  order  of  10-20  pm  [22],  If  can  be  shown  through  Mie  theory  [57],  given  the 
index  of  refraction,  number  density,  and  droplet  size,  that  these  particles  act  as  a  diffuse 
surface.  The  extinction  from  forward  scatfering  towards  the  head  face  prevenfs  reflecfed 
light  from  returning  to  the  camera,  and  the  angle  over  which  the  droplets  scatter  the 
incoming  light  is  large.  Therefore,  in  the  region  where  the  vapor  is  also  being  detected,  a 
strong  gradient  in  signal  is  present  where  the  liquid  to  vapor  transition  occurs.  The  rest 
of  the  cylinder  head,  which  is  more  diffusely  reflecfing  than  the  mirrored  section,  reflects 
much  more  of  the  incoming  LED  light  directly  back  to  the  camera  than  scattering  from  fhe 
droplef  cloud  forming  the  liquid  region.  Scattering  of  the  LED  light  by  the  dense  liquid 
droplets  attenuates  the  signal  from  its  typical  intensity  associated  with  reflection  off  the 
cylinder  head,  reducing  the  number  of  photons  returning  to  the  camera  in  regions  with 
liquid  fuel.  The  amount  of  backscatter  from  the  spray  toward  the  camera  is  less  than  the 
level  of  attenuation  caused  by  the  spray,  which  results  in  an  imaged  signal  with  lower 
intensity  than  the  surrounding  background  for  the  liquid  regions  in  the  spray. 

3.5  Post-Processing  Images 

As  discussed  in  the  literature  review,  much  of  the  recent  jet  measurement  work  uses 
single-pass  Schlieren  for  vapor  penetration  and  diffuse  back  illumination  (DBI)  for  liq¬ 
uid  penetration,  techniques  that  require  optical  access  from  both  sides  of  fhe  experiment 
[58].  Light  is  provided  from  one  side,  passes  through  the  jet,  and  is  imaged  on  the  other,  a 
setup  which  is  not  practical  in  an  optical  engine  when  simultaneously  measuring  liquid 
and  vapor  penetration  and  liftoff  length. 

An  effect  from  imaging  lighf  reflecfed  from  the  head  surface  is  that  the  background 
of  the  image  is  no  longer  uniform,  which  increases  the  complexity  in  distinguishing  the 
jet  intensity  from  the  background.  The  majority  of  the  literature  works  select  an  arbitrary 
pixel  intensity  threshold  that  best  bound  the  visible  jet  [8],  or  choose  a  threshold  centered 
between  peaks  in  an  intensity  histogram  [10],  methods  possible  for  images  with  constant 
backgrounds  where  the  jet  appears  as  a  uniformly  dark  region  against  a  bright  back¬ 
ground.  These  techniques  are  not  sufficient  for  a  non-consfanf  background,  as  shown  in 
parf  (a)  of  Figure  6,  where  the  jet  progresses  over  background  gradients  that  range  below 
and  above  the  maximum  intensity  values  of  the  jet. 

Images  taken  from  the  high-speed  camera  are  raster  images,  and  are  therefore  easily 
converted  to  three-dimensional  matrices  with  pixels  ranging  in  count  from  0  fo  2^,  where 
N  is  the  bit  depth  of  the  camera.  An  average  background  is  then  composed  from  a  mean 
of  the  ten  frames  preceding  the  frame  representing  SOI  and  subtracted  from  all  subse¬ 
quent  frames  to  be  analyzed.  The  target  frames  then  have  a  3x3  median  filter  applied 
to  them  to  reduce  shot  noise.  The  resultant  image  is  shown  in  Figure  6  (b).  The  images 
are  converted  to  polar  coordinates,  centered  on  the  injector  tip.  This  removes  the  need  to 
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(a)  (b)  (c)  (d) 

Figure  6:  Progressive  images  showing  points  in  the  post-processing  performed  on  high¬ 
speed  images  of  diesel  fuel  jets.  Frames  shown  from  (a)  original  video,  (b)  after  subtract¬ 
ing  background  and  applying  median  filter,  (c)  after  converting  to  polar  coordinates,  and 
(d)  converting  results  back  to  cartesian  coordinates. 
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track  the  jet  penetration  along  a  non-Cartesian  axis,  and  simplifies  the  dispersion  angle 
measurements. 


To  account  for  gradients  in  the  subtracted  non-constant  background  that  appear  in 
the  jet,  a  background  weighted  threshold  matrix  was  used  to  discern  between  true  back¬ 
ground  and  true  jet  signal.  First,  a  general  threshold  was  chosen  by  applying  Otsu's  statis¬ 
tical  method  to  the  final  background  subtracted  image,  where  the  jet  was  fully  developed. 
This  method  uses  a  closed-form  solution  to  minimize  the  average  error  in  assigning  pix¬ 
els  to  classes  of  average  intensity.  The  algorithm  maximizes  between  class  variance  by 
choosing  optimum  thresholds  for  each  class.  The  effectiveness  rj  of  this  method  is  calcu¬ 


lated  using  Equation  7: 


(7) 


where  us  is  the  between-class  variance  for  class  k,  and  ac  is  the  variance  of  the  entire 
image  [59].  The  number  of  classes  was  iteratively  increased  until  the  effectiveness  r]  ex¬ 
ceeded  99%.  The  general  threshold  was  then  set  as  the  mean  of  the  thresholds  associated 
with  the  two  lowest  values  of  k,  which  represent  the  edges  of  the  jet. 

A  matrix  the  size  of  the  polar  images  was  created,  filled  with  this  general  threshold 
value,  and  then  scaled  over  the  radial  axis  r  by  the  intensity  distribution  of  the  incoming 
power  of  the  LED  to  compensate  for  non-uniform  lighting.  The  polar  version  of  the  av¬ 
erage  background  was  then  subtracted  from  this  matrix,  resulting  in  a  threshold  matrix 
whose  level  is  chosen  statistically,  and  sensitive  to  non-constant  backgrounds  and  overall 
intensity  non-uniformity.  Edge  refinement  and  noise  rejection  code  is  applied  along  with 
the  threshold.  This  algorithm  has  performed  flexibly  over  different  lighting  and  back¬ 
ground  conditions  without  needing  to  tune  constants  or  modify  threshold  levels.  The 
liquid  was  distinguished  from  the  vapor  phase  in  a  similar  maimer,  by  differentiating  be¬ 
tween  the  Otsu  classes  in  the  two-phase  spatial  range  and  only  allowing  low  intensity 
classes  enclosed  by  vapor  to  be  classified  as  liquid. 

The  equation  defining  the  spray  angle,  9,  is  given  by  8: 


9  =  tan  ^ 


/  Ap,S/2  \ 

V(^/2)V 


(8) 


where  Ap^s/2  is  the  projected  spray  area  of  the  upstream  half  of  the  spray,  and  S  is  the 
penetration  length.  The  penetration  is  defined  as  the  distance  along  the  spray  axis  where 
1  /2  of  the  pixels  on  an  arc  of  6*/2  centered  on  the  axis  contain  signal  [10].  Determination  of 
the  regions  of  high-temperature  combustion  was  accomplished  using  the  Otsu  statistical 
technique  as  well,  although  the  videos  corresponding  to  the  OH  chemiluminescence  have 
a  uniform  background  at  the  level  of  the  dark  noise  of  the  camera,  removing  the  need  for 
a  spatially  weighted  threshold. 
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Figure  7:  Prcxress  of  aligning  penetration  data,  adjusting  for  variance  in  SOC-SOI  delay. 
Original,  uncorrected  data  sets  from  post-processing  algorithm  in  part  (a),  time-corrected 
sets  in  part  (b). 


3.6  Uncertainty  Quantilfication  and  Ensemble  Averaging 

The  uncertainty  of  single  measured  quantities  and  their  uncertainties  is  represented  by 
Equation  9: 

Xi  =  Aj  ±  AXi  (95%  confidence).  (9) 

The  value  of  Xi  is  described  by  a  measured  value  Xi^m  and  an  uncertainty  AXi  so  defined 
so  as  to  have  95%  confidence  of  Xi  falling  within  the  given  range.  Several  sets  of  opti¬ 
cal  data  were  taken  at  each  experimental  condition,  and  there  are  multiple  contributing 
sources  of  error  in  the  optical  measurements.  Therefore,  to  define  the  uncertainty  as  in¬ 
troduced  in  Equation  9  for  each  of  the  quantities  Xi  as  AXi,  another  root-sum-square  will 
be  used  to  combine  these  effects  [60]: 

AXi  =  ijAXli+AXl^+...  +  AXl,,.  (10) 

In  this  case,  AXi  represents  the  overall  uncertainty  for  the  given  quantity,  and 
represent  the  individual  uncertainty  components  making  up  the  total  uncertainty.  In  this 
work,  two  primary  sources  of  imcertainty  were  identified;  the  case-to-case  standard  de¬ 
viation  in  the  measurements  due  to  turbulence  effects  and  small  changes  in  experimental 
conditions,  and  uncertainty  due  to  the  experimental  setup  itself  due  to  the  limited  time 
and  spatial  resolutions  of  the  camera  and  optics  setup. .  The  uncertainty  interval  to  main¬ 
tain  95%  confidence  for  the  mean  (standard  deviation,  a  known)  of  a  multi-sample  experi¬ 
ment  where  the  true  mean  is  unknown  requires  interval  estimation  by  assuming  a  normal 
distribution  of  the  data  about  the  true  mean.  Confidence  intervals  for  many  experiments 
can  be  calculated  by  calculating  the  Z-statistic,  a  value  associated  with  the  area  under  the 
normal  distribution  curve.  For  experiments  with  <120  samples,  a  t-statistic  is  calculated 
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for  the  area  under  the  Student-t  distribution,  which  is  similar  to  the  normal  distribution, 
and  allows  for  the  estimation  of  the  true  mean  by  the  sample  mean.  The  measurements 
made  in  the  data  presented  in  this  report  involved  8-12  runs  at  each  condition,  so  the  t- 
statistic  was  used  for  confidence  intervals  [60,  61].  The  equation  for  the  t-statistic  is  given 
in  Equation  11: 


X-p 

s 

y/n 


(11) 


Here,  n  is  the  sample  size,  is  the  sample  variance  (square  of  sample  standard  deviation), 
X  is  the  sample  mean,  and  p  is  the  population  mean.  The  t-statistic  is  a  known  value  as 
a  function  of  %  confidence,  so  the  interval  about  a  population  mean  for  a  given  sample 
mean  and  standard  deviation  is  calculated  as: 


X  ± 

yn 


(12) 


The  data  presented  in  the  results  are  primarily  ensemble  averages  of  multiple  data 
sets.  However,  a  small  amount  of  fluctuation  in  the  time  between  electronic  start  of  injec¬ 
tion  command  (SOC)  and  the  actual  SOI  provided  an  incorrect  penetration  shape  during 
averaging,  as  noted  in  [25].  To  remove  the  uncertainty  associated  with  variability  in  the 
exact  SOI,  the  exact  timing  for  the  initial  ramp-up  of  the  penetration  rates  were  shifted  in 
time  such  that  the  root  mean  square  error  of  the  difference  between  the  initial  slopes  of 
the  ROI  of  each  case  was  minimized,  as  seen  in  Figure  7.  Data  presented  in  the  following 
section,  unless  otherwise  specified,  is  the  mean  of  injection  events  where  ranges 
between  8  and  10. 


4  Results 

Jet  penetration,  liquid  length,  jet  dispersion  angle,  liftoff  length,  and  engine  cylinder  pres¬ 
sure  data  have  been  acquired  over  a  range  of  compression-ignition  engine  operating  con¬ 
ditions  and  fuel  types.  In  this  section,  the  general  trends  of  these  data  over  the  range  of 
conditions  tested  are  first  presented  and  compared  to  results  seen  in  the  literature.  Next, 
behavior  of  these  parameters  is  investigated  during  the  initial  transient  development  pe¬ 
riod  of  the  jet  as  the  length  of  the  transient  period  is  varied.  Finally,  the  power  law  corre¬ 
lations  of  penetration  with  time  and  ambient  density  from  the  literature  are  compared  to 
the  correlations  drawn  to  represent  the  transient  data  from  the  current  work.  The  duration 
and  shape  of  the  unsteady  portion  of  the  ROI  was  a  function  of  the  injection  pressure  only. 
Plots  of  jet  penetration,  liquid  penetration,  and  dispersion  angle  are  bounded  in  time  by 
SOI  and  the  start  of  combustion.  Fiftotf  length  data  begins  when  high-temperature  signal 
was  first  detected,  and  ends  when  the  mean  liftoff  length  reached  the  injector  orifice  or 
too  much  signal  existed  in  the  image  to  correctly  capture  the  correct  length.  For  the  jet 
parameters,  data  markers  in  plots  are  visible  once  for  every  five  data  points,  and  once  for 
every  three  points  in  the  liftoff  plots. 
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4.1  Effects  of  Ambient  Density,  Injection  Pressure,  Intake  Tempera¬ 
ture,  and  Fuel  Type  on  Jet  Development 

High-speed  measurements  of  jet  and  combustion  parameters  were  acquired  for  more  than 
50  different  operating  conditions  encompassing  variations  in  injection  pressure,  fuel  type, 
intake  temperature,  and  ambient  cylinder  density  during  the  injection.  The  effects  of  these 
four  parameters  on  liftoff  length,  jet  penetration,  liquid  length,  and  dispersion  angle  are 
discussed,  and  example  figures  of  the  results  are  shown.  Eight  to  ten  runs  were  aver¬ 
aged  for  the  data  in  each  case  presented  in  the  data.  Error  bars  in  the  plots  represent  the 
standard  deviation  about  the  mean  of  the  samples  taken,  represented  by  markers.  The 
shaded  regions  represent  the  95%  confidence  interval  bands  for  the  mean  of  the  data. 
Start  of  combustion  times  as  found  through  analysis  of  the  mean  cylinder  pressure  traces 
are  shown  as  vertical  lines  on  the  liftoff  length  plots  and  cylinder  pressure/heat  release 
plots.  Methodology  for  finding  the  the  start  of  combustion  is  taken  from  [62].  In  direct- 
injection  compression  ignition  engines,  the  injection  of  liquid  fuel  into  the  hot  cylinder 
gases  causes  an  evaporative  cooling  effect,  resulting  in  a  drop  in  cylinder  pressure  rela¬ 
tive  to  an  equivalent  motored  trace.  The  ignition  delay  is  defined  as  the  location  where  the 
fired  cylinder  pressure  trace  recovers  its  motored  value  due  to  combustion  overcoming 
evaporative  cooling. 

4.1.1  Influence  of  Injection  Pressure 

The  length  and  shape  of  the  initial  transient  were  controlled  through  the  injection  pres¬ 
sure,  but  this  parameter  has  other  effects  on  jet  development  as  well.  Three  injection 
pressures  were  tested  at  each  set  of  engine  thermodynamic  conditions  and  fuels.  Time- 
resolved  results  for  the  selected  case  with  Ti^take  =  115°  C,  pamb  =  15  kg/m^  using  diesel 
fuel  are  shown  in  Eigure  8.  The  strong  correlation  between  injection  pressure  and  overall 
jet  penetration  is  evident.  As  the  fuel  pressure  increases,  the  increased  liquid  flow  veloci¬ 
ties  at  the  injector  orifice  translate  into  higher  penetration  rates. 

Conversely,  changes  in  injection  pressure  have  little  effect  on  the  quasi-steady  liquid 
length,  as  the  increase  in  injection  velocity  is  offset  by  the  increased  air  entrainment  rate 
[28].  However,  while  the  liquid  lengths  at  each  condition  reach  similar  mean  values,  the 
variance  in  the  data  decreases  with  increased  injection  pressure.  The  fluctuations  in  liq¬ 
uid  length  were  much  larger  for  lower  injection  pressures.  Eor  the  lower  injection  pressure 
cases,  isolated  regions  of  liquid  would  detach  from  the  end  of  spray,  causing  a  large  sud¬ 
den  drop  in  the  liquid  length,  a  phenomenon  less  prevalent  for  higher  fuel  pressures.  A 
sequence  of  images  illustrating  this  effect  is  shown  in  Eigure  10. 

Although  the  final  quasi-steady  value  of  the  dispersion  angle  for  all  three  injection 
pressures  was  similar,  the  lowest  injection  pressure  case  also  displayed  a  transient  dis¬ 
persion  angle  until  the  quasi-steady  period  of  the  rate-of-injection  profile  was  reached, 
approximately  375  ps  after  the  SOI.  This  increase  in  dispersion  angle  seems  to  correspond 
to  a  slight  decrease  in  the  liquid  length.  As  the  dispersion  angle  decreases  to  its  final 
quasi-steady  value  the  liquid  length  increases  to  its  quasi-steady  value.  This  effect  de¬ 
creased  with  increasing  injection  pressure,  to  a  point  where  the  dispersion  of  the  jet  at  the 
highest  injection  pressure  was  nearly  constant  throughout  the  entire  injection. 
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Figure  8:  Optical  results  for  a.)  jet  penetration,  b.)  liquid  penetration,  c.)  dispersion  angle, 
and  d.)  liftoff  length  against  time  for  all  injection  pressures  tested.  Intake  temperature  of 
115°C,  in-cylinder  density  of  15  kg/m^,  diesel  fuel. 


The  injection  pressure  has  a  strong  influence  on  the  ignition  delay  as  illustrated  by 
the  earlier  start  of  measurement  for  the  liftoff  length  as  injection  pressure  is  increased. 
Higher  injection  pressures  provided  shorter  ignition  delays.  These  trends  extend  to  the 
other  thermodynamic  cases  measured  as  well.  Although  not  as  obvious  as  the  influence 
of  injection  pressure  on  ignition  delay,  the  initial  liftoff  length  following  ignition  is  gener¬ 
ally  increased  with  higher  injection  pressure.  Following  ignition,  the  liftoff  length  for  all 
injection  pressures  rapidly  decreases  due  to  compressive  heating  of  the  in-cylinder  ambi¬ 
ent  gases.  The  ignition  delays  determined  from  the  cylinder  pressure  measurements  are 
consistently  2  CAD  (278  fis)  longer  than  the  delay  between  SOI  and  chemiluminescence 
initially  appearing  in  the  images.  The  chemiluminescence  first  observed  around  the  time 
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Figure  9:  Engine  cylinder  pressure  and  heat  release  rates  for  all  injection  pressures  tested. 
Intake  temperature  of  115°C,  in-cylinder  density  of  15  kg/  m^,  diesel  fuel. 
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Figure  10;  Sequence  of  images  showing  the  progression  of  the  liquid  length  and  subse¬ 
quent  detachment  of  a  Liquid  segment.  In  the  images,  the  vapor  boundary  is  outlined  in 
red  whereas  the  boundary  between  liquid  and  vapor  is  outlined  in  green. 


of  ignition  starts  as  small  kernels  that  grow  with  time  into  a  larger  region.  These  small 
kernels  are  trackable  optically  and  associated  with  the  liftoff  length,  but  it  is  not  until 
sufficient  heat  release  has  occurred  in-cylinder  to  overcome  the  pressure  decrease  due  to 
evaporative  cooling  that  the  ignition  delay  is  determined  from  the  in-cylinder  pressure 
measurement.  Tliis  results  in  the  observed  delay  between  the  chemiluminescence  based 
determination  of  ignition  delay  and  that  determined  from  the  in-cylinder  measured  pres¬ 
sure. 

4.1.2  Influence  of  Intake  Air  Temperature 

Jet  and  combustion  parameters  for  the  case  with  Pamb  =  15  kg/m^  and  Pi^j  =  30  MPa 
using  diesel  fuel  are  presented  in  Figure  11.  As  was  expected,  the  jet  penetration  was 
not  affected  significantly  by  changes  in  intake  temperature  when  the  ambient  density 
is  held  constant.  The  liquid  lengths  were  also  not  significantly  affected  by  this  change. 
Although  temperature  is  a  key  factor  in  the  vaporization  rate  of  the  fuel  jet,  the  range  of 
temperatures  swept  for  the  intake  conditions  was  not  large  enough  to  strongly  affect  the 
liquid  length  during  injection. 

This  range  of  temperatures  at  engine  intake  were  also  not  wide  enough  to  strongly 
affect  the  liftoff  lengths.  Liftoff  length  in  general  showed  much  larger  standard  deviation 
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Figure  11:  Optical  results  for  a.)  jet  penetration,  b.)  Liquid  penetration,  c.)  dispersion  angle, 
and  d.)  liftoff  length  against  time  for  aU  intake  temperatures  tested.  Injection  pressure  of 
30  MPa,  in-cylinder  density  of  15  kg/ m^,  diesel  fuel. 


than  the  other  parameters  measured.  This  variance  was  evident  in  both  the  location  and 
time  of  the  initial  ignition  kernel. 

4.1.3  Influence  of  In-Cylinder  Density 

Measurements  were  taken  for  two  in-cylinder  ambient  densities  specified  at  TDC,  which 
occurred  halfway  through  the  duration  of  the  injection.  In-cylinder  thermodynamic  con¬ 
ditions  in  an  reciprocating  engine  vary  with  time  due  to  piston  motion.  However,  the 
piston  speed  near  TDC  is  relatively  slow  and  the  injection  event  lasts  10  CAD,  from  -5 
to  +5  degrees  TDC,  causing  a  small  but  finite  density  change  of  3-4%  over  the  injection 
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Figure  12:  Engine  cylinder  pressure  and  heat  release  rates  for  all  intake  temperatures 
tested.  Injection  pressure  of  30  MPa,  in-cylinder  density  of  15  kg/ m^,  diesel  fuel. 
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event  for  all  tested  cases. 
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Figure  13:  Optical  results  for  a.)  jet  penetration,  b.)  liquid  penetration,  c.)  dispersion  angle, 
and  d.)  liftoff  length  against  time  for  all  in-cylinder  densities  tested,  and  diesel  fuel.  Intake 
temperature  of  115°C,  injection  pressure  of  30  MPa. 


Ambient  density  was  seen  to  reduce  the  quasi-steady  values  of  liquid  length  and  va¬ 
por  penetration  length,  as  seen  in  the  case  using  diesel  fuel  with  an  intake  temperature 
of  115°C,  and  injection  pressure  of  30  MPa  shown  in  Figure  13.  This  is  as  expected,  jet 
penetration  and  liquid  length  scale  inversely  in  a  non-linear  manner  with  ambient  den¬ 
sity  in  the  long  time  limit  [11,  28].  For  short  times,  the  jet  penetration  was  slightly  larger 
for  higher  ambient  density  cases,  as  seen  in  Figures  13  and  14.  This  effect  was  greater  for 
diesel  fuel,  and  should  be  explored  further.  The  dispersion  angle  did  not  show  any  statis¬ 
tically  significant  trends  with  density  once  the  quasi-steady  injection  period  was  reached. 

Liftoff  lengths  were  strongly  affected  by  changes  in  ambient  density.  The  overall  mean 
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Figure  14:  Additonal  jet  penetration  results  over  all  densities  for  a.)  JP8  and  b.)  HRJ.  Intake 
temperature  of  115°C,  injection  pressure  of  30  MPa. 


liftoff  length  decreased  with  increasing  density,  and  the  ignition  delay  was  severely  re¬ 
duced  with  the  increase  in  ambient  density.  Increasing  ambient  density  enhances  jet  mix¬ 
ing,  resulting  in  faster  energy  transfer  from  hot  cylinder  gases  to  the  fuel. 

4.1.4  Influence  of  Fuel  Type 

The  jet  fuels  had  higher  volatility  than  the  #2  diesel  fuel,  decreasing  the  distance  and 
time  required  for  the  liquid  fuel  to  evaporate  completely.  This  is  reflected  in  data  for  case 
Pamb  =  15  kg/m^  and  Pi„j  =  30  MPa  in  Figure  16.  The  jet  fuels  displayed  shorter  jet 
penetration  and  liquid  lengths  than  the  diesel  test  fuel  due  to  the  volatility  differences 
between  the  fuels. 

The  combustion  results  did  not  follow  the  same  trend.  The  diesel  and  JP8  fuels  de¬ 
veloped  similar  liftoff  and  ignition  delay  profiles,  while  the  HRJ  fuel  ignited  earlier  and 
maintained  a  liftoff  distance  closer  to  the  injection  orifice.  As  shown  in  Table  2,  HRJ  has  a 
significantly  higher  cetane  number  than  the  other  two  fuels  indicating  it  should  have  the 
shortest  ignition  delay,  a  result  seen  in  Figure  17.  The  #2  diesel  fuel  and  the  JP8  fuel  used 
have  cetane  numbers  that  are  the  same  within  the  reproducibility  limits  of  the  cetane 
number  test.  The  matched  liftoff  length  history  and  ignition  delays  for  the  JP8  and  #2 
diesel  fuel  agree  with  previous  observations  that  the  volatility  difference  between  the  fu¬ 
els  does  not  significant  impact  ignition  delay  and  combustion  if  cetane  number  is  matched 
between  the  fuels  [62,  63].  The  results  are  reflected  in  the  ensemble  averaged  engine  cylin¬ 
der  pressure  results  and  associated  heat  release  results  presented  in  Figure  17.  The  heat 
release  and  cylinder  pressure  results  confirm  the  similarity  of  the  combustion  process  for 
the  #2  diesel  and  JP8  fuels.  It  is  evident  from  the  heat  release  that  for  all  fuels,  a  large 
percentage  of  the  fuel  is  combusted  in  the  initial  premixed  bum,  with  a  relatively  small 
fraction  of  heat  release  associated  with  the  mixing  controlled  combustion  at  these  low 
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Figure  15:  Engine  cylinder  pressure  and  heat  release  rates  for  all  in-cylinder  densities 
tested.  Intake  temperature  of  115°C,  injection  pressure  of  30  MPa,  diesel  fuel. 
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Figure  16:  Optical  results  for  a.)  jet  penetration,  b.)  liquid  penetration,  c.)  dispersion  tingle, 
and  d.)  liftoff  length  against  time  for  all  fuels  tested.  Intake  temperature  of  115°C,  injection 
pressure  of  30  MPa,  in-cylinder  density  of  15  kg/m^. 


density  conditions  where  the  ignition  delay  is  long. 

4.2  Impact  of  Initial  Transient  ROI  on  Jet  Development 

It  was  shown  in  the  previous  section  that  in  an  optically  accessible  internal  combustion 
engine  ,  during  the  quasi-steady  portion  of  the  injection  process  the  jet  and  combustion 
parameters  measured  agreed  with  general  trends  from  similar  experiments  in  constant 
volume  chambers  [10,  28].  The  following  section  explores  deviations  from  quasi-steady 
behavior  when  the  injection  velocity  is  not  constant  during  the  initial  ramp-up  of  the 
injection  process. 
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Figure  17:  Engine  cylinder  pressure  and  heat  release  rates  for  all  fuels  tested.  Intake  tem¬ 
perature  of  115°C,  injection  pressure  of  30  MPa,  in-cylinder  density  of  15  kg/m^. 
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4.2.1  Jet  Behavior  During  Transient  ROI 

The  transient  part  of  the  injection  event  is  due  to  a  varying  fuel  flow  rate  as  the  injec¬ 
tor  needle  lifts  from  covering  the  injector  orifices.  As  the  fuel  pressure  overcomes  the 
spring  force  holding  the  injector  shut,  the  flow  area  increases  to  its  maximum  value  over 
a  non-negligible  time  duration,  and  the  mass  flow  rate  changes  from  initially  zero  to  a 
quasi-steady  value.  This  rise  in  mass  flow  rate  affects  the  penetration  rate  of  the  jet  dur¬ 
ing  its  initial  development  when  it  is  dominated  by  the  momentum  of  the  liquid  spray. 
The  penetration  relationship  proposed  by  Hiroyasu  [11]  predicts  a  linear  scaling  of  pene¬ 
tration  with  time  until  hreak,  then  a  scaling  afterwards.  For  this  work,  the  best  tit  for 
the  scaling  of  jet  penetration  with  time  was  found.  Similar  to  previous  work,  a  power  law 
equation  of  the  form: 

S{t)  =  (13) 

was  found  to  fit  the  data  well,  where  /I  is  a  empirical  constant  and  a  is  the  power  law 
exponent.  Three  values  of  a  and  /3  were  found  for  each  condition  tested,  one  pair  fit  to  the 
initial  (short  time)  data  points,  one  to  the  transition  period  of  the  jet,  and  the  other  tit  from 
data  in  the  quasi-steady  (long-time)  regime.  The  jet  penetration  data  initially  presented  in 
Figure  8  is  shown  on  a  log-log  scale,  overlaid  with  the  curve  fits  given  by  Equation  13,  is 
in  Figure  18.  The  values  of  a  derived  to  best  tit  the  data  are  also  included  in  the  plot. 

Three  linear  segments  present  themselves  in  jet  penetration  on  a  log-log  scale,  sepa¬ 
rated  by  clear  inflection  points.  Similar  to  previous  results  in  the  literature,  the  depen¬ 
dence  on  time,  a,  for  the  long-time  limits  were  all  close  to  a  p/2 

scaling.  However,  in 

the  period  before  the  jet  transitions  from  its  ramp-up  state  while  the  opening  transient  is 
present,  the  experimental  results  deviate  from  the  literature.  Equation  1  predicts  a  linear 
relationship  with  time  during  this  period,  but  the  fits  derived  from  the  data  displayed 
a  nonlinear  scaling  with  time.  This  results  was  seen  for  measurements  at  all  conditions 
tested.  The  dependence  from  the  literature  assumes  a  constant  mass  flow  velocity,  and 
during  the  liquid-momentum  dominant  phase,  the  penetration  velocity  would  also  be 
constant.  The  acceleration  in  mass  flow  rate  due  to  the  opening  of  the  injector  tip,  how¬ 
ever,  increases  this  dependence  to  be  >1,  as  seen  in  Figure  18. 

It  is  clear  that  S  is  strongly  affected  by  the  opening  of  the  injector  during  the  ini¬ 
tial  development  stage,  but  transitions  to  conventional  quasi-steady  behavior  afterwards. 
Therefore,  although  the  initial  ROI  transient  affects  the  jet  penetration  while  the  transient 
is  occurring,  the  head  of  the  jet  eventually  transitions  to  the  time  dependence  predicted  in 
the  literature  [11].  Although  quasi-steady  dependence  of  penetration  on  time  is  achieved 
for  long  injection  durations,  a  significant  error  will  still  exist  due  to  the  difference  at  early 
times. 

4.2.2  Jet  Velocity  and  Unsteady  ROI  During  the  Injection  Transient 

In  order  to  directly  compare  the  effect  of  the  rate  of  fuel  injection  on  jet  penetration,  the 
velocity  of  the  jet  tip,  is  presented  with  the  time  rate  of  change  of  mass  flow,  in 
Figure  19.  The  derivative  of  the  curve  fit  defined  in  Equation  13  is  defined  by  the  following 
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Figure  18:  Jet  penetration  on  a  log-log  scale,  with  curve  fits  of  the  form  in  Equation  13 
included.  Diesel  fuel  at  an  intake  temperature  of  115°C,  and  in-cylinder  density  of  15 
kg/m^. 
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equation: 


(14) 


Figure  19:  Comparison  of  the  jet  penetration  velocity,  rate  of  ROI,  and  the  derivative  of 
the  curve  fit  from  Equation  13.  Conditions  of  Pamh-^^  kg/m^,  T'i„tafce=115°C,  using  diesel 
fuel. 

The  jet  penetration  velocity  as  predicted  by  Equation  14  correctly  represents  the  bound¬ 
ary  condition  of  ^  0  for  f  ^  0.  This  is  in  contrast  with  the  non-physical  constant  spray 

tip  velocity  down  to  t=0  implied  from  the  derivative  of  Equation  1,  an  issue  noted  in  [25]. 
Excluding  the  period  of  transition  where  the  penetration  shifted  from  its  short  to  long 
time  behavior,  the  derivative  of  the  curve  fit  describes  the  behavior  of  the  actual  spray  tip 
velocity  reasonably  well. 

The  shape  of  also  appears  to  follow  a  similar  trend  as  the  time  dependence  of  the 
jet  velocity  during  the  transient  period.  The  correlation  between  ROI  and  tip  penetration 
velocity  falls  off  as  the  jet  reaches  its  long-time  behavior.  This  relationship  between  the 
derivative  of  the  ROI  and  the  jet  tip  penetration  velocity  provides  an  possible  method  of 
approximating  jet  penetration  for  cases  with  injection  transients.  Rather  than  assuming  a 
constant  mass  flow  rate,  the  actual  time  dependence  of  the  ROI  of  the  injector  could  be 
measured  and  possible  relations  between  initial  ROI  and  penetration  developed.  Injector 
ROI  is  a  much  easier  quantity  to  experimentally  measure  than  actual  optical  jet  data. 

4.2.3  Transient  ROI  Effects  on  Liftoff  Length  and  Combustion 

The  experimental  conditions  run  were  such  that  the  delay  between  SOI  and  initial  com¬ 
bustion  were  relatively  long  compared  to  the  time  between  SOI  and  the  shift  of  the  jet  to 
its  quasi-steady  region.  Therefore,  by  the  time  high-temperature  reactions  began,  the  jet 
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was  no  longer  experiencing  the  effects  of  the  initial  transient.  The  main  effect  the  non¬ 
constant  ROI  has  on  the  combustion  of  jets  under  these  conditions  is  the  change  in  fuel 
mass  injecfed.  If  the  injection  event  truly  followed  an  instant-on  profile,  the  mass  injected 
would  be  greater,  resulting  in  a  change  in  engine  load. 

Transients  of  this  form,  then,  exert  secondary  effects  on  the  successive  combustion 
event.  For  instance,  the  transient  injection  will  alter  the  rate  of  air  entrainment  and,  there¬ 
fore,  the  subsequent  ignition  delay  and  initially  established  liftoff  length.  Again,  the  ac¬ 
tual  effect  the  transient  has  on  combustion  is  increased  as  the  injection  duration  decreases. 
For  cases  where  the  injection  ends  before  it  reaches  its  quasi-steady  state,  the  effects  of  the 
transient  on  fuel  mass  and  mixing  will  more  directly  affect  the  ignition  location  and  liftoff 
length. 


5  Summary  and  Impact 

As  mentioned  in  the  results,  in  the  environment  of  an  internal  combustion  engine,  the 
quasi-steady  regions  of  the  jet  and  combustion  parameters  measured  agreed  with  general 
trends  seen  in  scaling  predictions  from  constant  volume  chambers  in  the  literature.  This 
indicates  that  with  the  correct  tuning  parameters,  the  relationships  developed  between 
the  testing  conditions  and  the  resulting  jet  features  for  chambers  can  be  applied  in  en¬ 
gine  work  as  well.  However,  it  is  now  evident  that  this  is  true  only  for  cases  where  the 
the  injection  is  dominated  by  its  quasi-steady  region.  This  happens  when  either  the  the 
transient  time  is  short  relative  to  the  steady  portion,  or  if  the  injection  follows  an  instanf- 
on,  "top-hat"  style  of  ramp-on.  The  error  caused  by  deviations  from  these  scaling  laws 
increases  as  the  ratio  of  transient  to  non- transient  portions  increases. 

The  primary  results  and  conclusions  based  on  the  findings  from  the  current  study  are: 

•  High-speed  measurements  of  jet  parameters  and  compression  ignition  combustion 
were  obtained  simultaneously  in  an  optical  engine  whose  specifications  are  similar 
to  those  of  modem  commercial  small-bore  diesel  engines.  This  was  achieved  by  de¬ 
signing  and  implementing  a  front-illuminated  optical  setup  that  used  Mie  scattering 
for  liquid  length  determination  and  dispersion  angle  measurements,  double-pass 
shadow graphy  for  vapor  measurements,  and  chemiluminescence  for  determination 
of  liftoff  length.  The  jet  measurements  were  obtained  at  a  frame  rate  of  121  kHz  so 
as  to  capture  in  detail  the  penetration  during  the  entire  injection,  with  key  interest 
in  the  initial  ramp-up  transient  seen  in  hydraulically-activated  fuel  injectors. 

•  High-pressure  fuel  jets  in  internal  combustion  engines  were  found  to  follow  the 
same  trends  with  ambient  and  injection  conditions  as  results  taken  previously  in 
constant  volume  chambers.  Key  observations  include  reduction  in  jet  and  liquid 
penetration  with  ambient  density,  increase  in  jet  penetration  with  injection  pressure, 
and  decrease  in  quasi-steady  liquid  length  with  fuel  volatility. 

•  The  influence  of  fuel  properties  on  the  jet  parameters  was  also  characterized  using 
two  jet  fuels.  The  higher  volatility  and  lower  density  of  the  jet  fuels  relative  to  #2 
diesel  fuel  resulted  in  shorter  liquid  lengths  for  the  jet  fuels  and  slight  differences  in 
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the  jet  dispersion  and  penetration.  However,  the  lift-off  length  and  ignition  delay  for 
the  JP8  fuel  were  almost  identical  to  the  #2  diesel  fuel  with  matched  cetane  number, 
indicating  that  cetane  number  has  a  stronger  influence  than  changes  in  fuel  physical 
properties  on  the  combustion  process. 

•  The  effects  of  the  injector  opening  process  were  studied  by  varying  the  injection 
pressure,  which  modified  the  fuel  mass  rate-of-injection  from  initially  zero  to  a 
quasi-steady  value.  Analysis  of  the  jet  penetration  results  revealed  three  distinct 
scalings  with  time,  as  opposed  to  the  two-regime  methodology  used  in  the  litera¬ 
ture.  The  jet  was  seen  to  initially  have  a  scaling  with  time  of  with  a  >  1,  moving 
to  a  stage  where  the  scaling  was  approximately  linear,  and  eventually  transitioning 
to  the  scaling  expected  for  these  types  of  jets  in  the  long-time  domain  [11].  The 
initial  time  scaling  varied  with  injection  pressure,  and  the  transition  to  ^1/2 

scaling 

with  time  indicates  the  point  at  which  the  transient  has  ended,  and  the  jet  is  follow¬ 
ing  the  scaling  predicted  for  long  times  from  the  literature  [11]. 

•  The  penetration  velocity  of  the  jet  was  compared  to  the  derivative  of  the  ROI,  or  the 
rate  of  change  of  mass  flow  from  the  injector.  For  the  cases  considered,  these  param¬ 
eters  followed  similar  scaling  with  time.  This  suggests  that  the  ROI  of  a  jet,  which  is 
experimentally  much  easier  to  measure  than  actual  optical  penetration  data,  could 
possibly  inform  transient  effects  on  the  scaling  of  jet  parameters  with  time,  although 
further  work  is  needed  to  confirm  this  theory. 

The  transients  introduced  by  the  injector  opening  process  influenced  the  jet  tip  pene¬ 
tration  velocity  for  the  duration  of  the  opening  process.  Jet  development  then  transitioned 
to  behavior  expected  of  a  quasi-steady  process  after  the  injector  was  fully  open.  There¬ 
fore,  injections  whose  quasi-steady  portion  is  long  relative  to  the  transient  period  can 
be  well  described  by  the  scaling  relationships  from  the  literature,  while  short  injections 
dominated  primarily  by  the  opening  and  closing  periods  of  the  injector  cannot  assume 
an  instant-on  injection  profile,  and  must  take  the  rate  of  injection  scaling  with  time  into 
consideration  to  correctly  model  jet  behavior. 
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